Estuarine frontogenesis

Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 546–561, doi:10.1175/JPO-D-14-0082.1.Model studies and observations in the Hudson River estuary indicate that frontogenesis occurs as a result of topographic forcing. Bottom fronts form just downstream of lateral constrictions, where the width of the estuary increases in the down-estuary (i.e., seaward) direction. The front forms during the last several hours of the ebb, when the combination of adverse pressure gradient in the expansion and baroclinicity cause a stagnation of near-bottom velocity. Frontogenesis is observed in two dynamical regimes: one in which the front develops at a transition from subcritical to supercritical flow and the other in which the flow is everywhere supercritical. The supercritical front formation appears to be associated with lateral flow separation. Both types of fronts are three-dimensional, with strong lateral gradients along the flanks of the channel. During spring tide conditions, the fronts dissipate during the flood, whereas during neap tides the fronts are advected landward during the flood. The zone of enhanced density gradient initiates frontogenesis at multiple constrictions along the estuary as it propagates landward more than 60 km during several days of neap tides. Frontogenesis and frontal propagation may thus be essential elements of the spring-to-neap transition to stratified conditions in partially mixed estuaries.Support for this research was provided by NSF Grant OCE 0926427.2015-08-0

Sediment transport time scales and trapping efficiency in a tidal river

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Earth Surface 122 (2017): 2042–2063, doi:10.1002/2017JF004337.Observations and a numerical model are used to characterize sediment transport in the tidal Hudson River. A sediment budget over 11 years including major discharge events indicates the tidal fresh region traps about 40% of the sediment input from the watershed. Sediment input scales with the river discharge cubed, while seaward transport in the tidal river scales linearly, so the tidal river accumulates sediment during the highest discharge events. Sediment pulses associated with discharge events dissipate moving seaward and lag the advection speed of the river by a factor of 1.5 to 3. Idealized model simulations with a range of discharge and settling velocity were used to evaluate the trapping efficiency, transport rate, and mean age of sediment input from the watershed. The seaward transport of suspended sediment scales linearly with discharge but lags the river velocity by a factor that is linear with settling velocity. The lag factor is 30–40 times the settling velocity (mm s−1), so transport speeds vary by orders of magnitude from clay (0.01 mm s−1) to coarse silt (1 mm s−1). Deposition along the tidal river depends strongly on settling velocity, and a simple advection-reaction equation represents the loss due to settling on depositional shoals. The long-term discharge record is used to represent statistically the distribution of transport times, and time scales for settling velocities of 0.1 mm s−1 and 1 mm s−1 range from several months to several years for transport through the tidal river and several years to several decades through the estuary.Hudson River Foundation Grant Number: 004/13A; National Science Foundation Grant Number: 13251362018-05-0

Response to channel deepening of the salinity intrusion, estuarine circulation, and stratification in an urbanized estuary

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research-Oceans 124(7), (2019): 4784-4802, doi: 10.1029/2019JC015006.Modifications for navigation since the late 1800s have increased channel depth (H) in the lower Hudson River estuary by 10–30%, and at the mouth the depth has more than doubled. Observations along the lower estuary show that both salinity and stratification have increased over the past century. Model results comparing predredging bathymetry from the 1860s with modern conditions indicate an increase in the salinity intrusion of about 30%, which is roughly consistent with the H5/3 scaling expected from theory for salt flux dominated by steady exchange. While modifications including a recent deepening project have been concentrated near the mouth, the changes increase salinity and threaten drinking water supplies more than 100 km landward. The deepening has not changed the responses to river discharge (Qr) of the salinity intrusion (~Qr−1/3) or mean stratification (Qr2/3). Surprisingly, the increase in salinity intrusion with channel deepening results in almost no change in the estuarine circulation. This contrasts sharply with local scaling based on local dynamics of an H2 dependence, but it is consistent with a steady state salt balance that allows scaling of the estuarine circulation based on external forcing factors and is independent of depth. In contrast, the observed and modeled increases in stratification are opposite of expectations from the steady state balance, which could be due to reduction in mixing with loss of shallow subtidal regions. Overall, the mean shift in estuarine parameter space due to channel deepening has been modest compared with the monthly‐to‐seasonal variability due to tides and river discharge.Funding was provided by NSF Coastal SEES (OCE 1325136). Data supporting this study are posted to Zenodo (https://doi.org/10.5281/zenodo.2551285) or are available by contacting the author.2019-12-0

Shear and turbulence production across subtidal channels

Author Posting. © Sears Foundation for Marine Research, 2006. This article is posted here by permission of Sears Foundation for Marine Research for personal use, not for redistribution. The definitive version was published in Journal of Marine Research 64 (2006): 147-171, doi:10.1357/002224006776412359.In intertidal regions with subtidal channels, effects of bathymetry on overlying flow vary greatly with tidal stage. Around low water when mudflats and marsh are exposed, flow is constrained to channels, but when water depths are greater, tidal forcing may not necessarily be aligned with meandering channel axes. Flow across the channel can generate strong shear and turbulence at the elevation of the channel banks and can significantly increase turbulent energy in the middle of the water column. Field observations in a mudflat channel of San Francisco Bay indicate that cross-channel shear regularly occurs there early in ebb tides. With increased freshwater flow, baroclinic forcing can enhance shear by decoupling flow between dense water flooding in the channel and fresher water ebbing above the channel banks. A water column numerical model with κ-ε turbulence closure is modified to represent the cross-channel shear production. Numerical results with uniform density indicate that turbulence production increases with the angle between the barotropic tidal forcing and the channel axis. When a longitudinal salinity gradient is imposed, cross-channel shear production contributes to breakdown of periodic stratification. Turbulence produced at the channel banks locally exceeds dissipation, and the excess energy is either lost to buoyancy or diffuses vertically to lower energy regions near the surface and near the bed. The balance among shear production, buoyancy production, and diffusion of turbulence depends on the flow angle and the strength of the longitudinal salinity gradient.This research was funded by National Institutes of Health grant no. P42ES0475 from the National Institute of Environmental Health Sciences

Tidal and meteorological forcing of sediment transport in tributary mudflat channels

Author Posting. © Elsevier B.V., 2007. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Continental Shelf Research 27 (2007): 1510-1527, doi:10.1016/j.csr.2007.01.010.Field observations of flow and sediment transport in a tributary channel through intertidal mudflats indicate that suspended sediment was closely linked to advection and dispersion of a tidal salinity front. During calm weather when tidal forcing was dominant, high concentrations of suspended sediment advected up the mudflat channel in the narrow region between salty water from San Francisco Bay and much fresher runoff from the small local watershed. Salinity and suspended sediment dispersed at similar rates through each tidal inundation, such that during receding ebbs the sediment pulse had spread spatially and maximum concentrations had decreased. Net sediment transport was moderately onshore during the calm weather, as asymmetries in stratification due to tidal straining of the salinity front enhanced deposition, particularly during weaker neap tidal forcing. Sediment transport by tidal forcing was periodically altered by winter storms. During storms, strong winds from the south generated wind waves and temporarily increased suspended sediment concentrations. Increased discharge down the tributary channels due to precipitation had more lasting impact on sediment transport, supplying both buoyancy and fine sediment to the system. Net sediment transport depended on the balance between calm weather tidal forcing and perturbations by episodic storms. Net transport in the tributary channel was generally off-shore during storms and during calm weather spring tides, and on-shore during calm weather neap tides.The research was funded by National Institutes of Health grant P42ES0475 from the National Institute of Environmental Health Sciences. Reference

Estuarine exchange flow variability in a seasonal, segmented estuary

Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(3),(2020): 595-613, doi:10.1175/JPO-D-19-0108.1.Small estuaries in Mediterranean climates display pronounced salinity variability at seasonal and event time scales. Here, we use a hydrodynamic model of the Coos Estuary, Oregon, to examine the seasonal variability of the salinity dynamics and estuarine exchange flow. The exchange flow is primarily driven by tidal processes, varying with the spring–neap cycle rather than discharge or the salinity gradient. The salinity distribution is rarely in equilibrium with discharge conditions because during the wet season the response time scale is longer than discharge events, while during low flow it is longer than the entire dry season. Consequently, the salt field is rarely fully adjusted to the forcing and common power-law relations between the salinity intrusion and discharge do not apply. Further complicating the salinity dynamics is the estuarine geometry that consists of multiple branching channel segments with distinct freshwater sources. These channel segments act as subestuaries that import both higher- and lower-salinity water and export intermediate salinities. Throughout the estuary, tidal dispersion scales with tidal velocity squared, and likely includes jet–sink flow at the mouth, lateral shear dispersion, and tidal trapping in branching channel segments inside the estuary. While the estuarine inflow is strongly correlated with tidal amplitude, the outflow, stratification, and total mixing in the estuary are dependent on the seasonal variation in river discharge, which is similar to estuaries that are dominated by subtidal exchange flow.We thank two anonymous reviewers for constructive comments, the staff of the South Slough National Estuarine Research Reserve for providing time series data, and Parker MacCready for sharing LiveOcean boundary conditions. This work was partially sponsored by the National Estuarine Research Reserve System Science Collaborative, which supports collaborative research that addresses coastal management problems important to the reserves. The Science Collaborative is funded by the National Oceanic and Atmospheric Administration and managed by the University of Michigan Water Center (NAI4NOS4190145). Computations were performed on the University of Oregon high performance computer Talapas.2020-08-2

Watershed Suspended Sediment Supply and Potential Impacts of Dam Removals for an Estuary

Observations and modeling are used to assess potential impacts of sediment releases due to dam removals on the Hudson River estuary. Watershed sediment loads are calculated based on sediment-discharge rating curves for gauges covering 80% of the watershed area. The annual average sediment load to the estuary is 1.2 Mt, of which about 0.6 Mt comes from side tributaries. Sediment yield varies inversely with watershed area, with regional trends that are consistent with substrate erodibility. Geophysical and sedimentological surveys in seven subwatersheds of the Lower Hudson were conducted to estimate the mass and composition of sediment trapped behind dams. Impoundments were classified as (1) active sediment traps, (2) run-of-river sites not actively trapping sediment, and (3) dammed natural lakes and spring-fed ponds. Based on this categorization and impoundment attributes from a dam inventory database, the total mass of impounded sediment in the Lower Hudson watershed is estimated as 4.9 ± 1.9 Mt. This represents about 4 years of annual watershed supply, which is small compared with some individual dam removals and is not practically available given current dam removal rates. More than half of dams impound drainage areas less than 1 km2, and play little role in downstream sediment supply. In modeling of a simulated dam removal, suspended sediment in the estuary increases modestly near the source during discharge events, but otherwise effects on suspended sediment are minimal. Fine-grained sediment deposits broadly along the estuary and coarser sediment deposits near the source, with transport distance inversely related to settling velocity

Waves in the Red Sea : response to monsoonal and mountain gap winds

Author Posting. © The Author(s), 2013. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Continental Shelf Research 65 (2013): 1-13, doi:10.1016/j.csr.2013.05.017.An unstructured grid, phase-averaged wave model forced with winds from a high resolution atmospheric model is used to evaluate wind wave conditions in the Red Sea over an approximately 2-year period. The Red Sea lies in a narrow rift valley, and the steep topography surrounding the basin steers the dominant wind patterns and consequently the wave climate. At large scales, the model results indicated that the primary seasonal variability in waves was due to the monsoonal wind reversal. During the winter, monsoon winds from the southeast generated waves with mean significant wave heights in excess of 2 m and mean periods of 8 s in the southern Red Sea, while in the northern part of the basin waves were smaller, shorter period, and from northwest. The zone of convergence of winds and waves typically occurred around 19-20˚N, but the location varied between 15 to 21.5˚N. During the summer, waves were generally smaller and from the northwest over most of the basin. While the seasonal winds oriented along the axis of the Red Sea drove much of the variability in the waves, the maximum wave heights in the simulations were not due to the monsoonal winds but instead were generated by localized mountain wind jets oriented across the basin (roughly east-west). During the summer, a mountain wind jet from the Tokar Gap enhanced the waves in the region of 18 and 20˚N, with monthly mean wave heights exceeding 2 m and maximum wave heights of 14 m during a period when the rest of the Red Sea was relatively calm. Smaller mountain gap wind jets along the northeast coast created large waves during the fall and winter, with a series of jets providing a dominant source of wave energy during these periods. Evaluation of the wave model results against observations from a buoy and satellites found that the spatial resolution of the wind model significantly affected the quality of the wave model results. Wind forcing from a 10-km grid produced higher skills for waves than winds from a 30-km grid, largely due to under-prediction of the mean wind speed and wave height with the coarser grid. The 30-km grid did not resolve the mountain gap wind jets, and thus predicted lower wave heights in the central Red Sea during the summer and along the northeast coast in the winter.This research is based on work supported by Award No. USA00001, USA00002, KSA00011, made by the King Abdullah University of Science and Technology (KAUST) in Saudi Arabia

Hydraulics and mixing in a laterally divergent channel of a highly stratified estuary

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 122 (2017): 4743–4760, doi:10.1002/2016JC012455.Estuarine mixing is often intensified in regions where topographic forcing leads to hydraulic transitions. Observations in the salt-wedge estuary of the Connecticut River indicate that intense mixing occurs during the ebb tide in regions of supercritical flow that is accelerated by lateral expansion of the channel. The zones of mixing are readily identifiable based on echo-sounding images of large-amplitude shear instabilities. The gradient Richardson number (Ri) averaged across the mixing layer decreases to a value very close to 0.25 during most of the active mixing phase. The along-estuary variation in internal Froude number and interface elevation are roughly consistent with a steady, inviscid, two-layer hydraulic representation, and the fit is improved when a parameterization for interfacial stress is included. The analysis indicates that the mixing results from lateral straining of the shear layer, and that the rapid development of instabilities maintains the overall flow near the mixing threshold value of Ri = 0.25, even with continuous, active mixing. The entrainment coefficient can be estimated from salt conservation within the interfacial layer, based on the finding that the mixing maintains Ri = 0.25. This approach leads to a scaling estimate for the interfacial mixing coefficient based on the lateral spreading rate and the aspect ratio of the flow, yielding estimates of turbulent dissipation within the pycnocline that are consistent with estimates based on turbulence-resolving measurements.NSF Grant Number: OCE 0926427; Devonshire Scholars program2017-12-1

Stratified turbulence and mixing efficiency in a salt wedge estuary

Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 1769-1783, doi:10.1175/JPO-D-15-0193.1.High-resolution observations of velocity, salinity, and turbulence quantities were collected in a salt wedge estuary to quantify the efficiency of stratified mixing in a high-energy environment. During the ebb tide, a midwater column layer of strong shear and stratification developed, exhibiting near-critical gradient Richardson numbers and turbulent kinetic energy (TKE) dissipation rates greater than 10−4 m2 s−3, based on inertial subrange spectra. Collocated estimates of scalar variance dissipation from microconductivity sensors were used to estimate buoyancy flux and the flux Richardson number Rif. The majority of the samples were outside the boundary layer, based on the ratio of Ozmidov and boundary length scales, and had a mean Rif = 0.23 ± 0.01 (dissipation flux coefficient Γ = 0.30 ± 0.02) and a median gradient Richardson number Rig = 0.25. The boundary-influenced subset of the data had decreased efficiency, with Rif = 0.17 ± 0.02 (Γ = 0.20 ± 0.03) and median Rig = 0.16. The relationship between Rif and Rig was consistent with a turbulent Prandtl number of 1. Acoustic backscatter imagery revealed coherent braids in the mixing layer during the early ebb and a transition to more homogeneous turbulence in the midebb. A temporal trend in efficiency was also visible, with higher efficiency in the early ebb and lower efficiency in the late ebb when the bottom boundary layer had greater influence on the flow. These findings show that mixing efficiency of turbulence in a continuously forced, energetic, free shear layer can be significantly greater than the broadly cited upper bound from Osborn of 0.15–0.17.Holleman was supported by the Devonshire Scholars program. The field study and the coauthors’ contributions were supported by NSF Grant OCE 0926427.2016-11-2